Statistics of Wave Functions in Coupled Chaotic Systems 
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Using the supersymmetry technique, we calculate the joint 
distribution of local densities of electron wavefunctions in two 
coupled disordered or chaotic quantum billiards. We find 
novel spatial correlations that are absent in a single chaotic 
system. Our exact result can be interpreted for small coupling 
in terms of the hybridization of eigenstates of the isolated 
billiards. We show that the presented picture is universal, 
independent of microscopic details of the coupling. 

PACS numbers: 05.45.Mt, 73.23.-b, 73.20.Dx 



In disordered or chaotic systems the dynamics is in 
general too complex to find exact eigenstates for a given 
disorder or boundary configuration. Therefore, one is 
interested in statistical properties of physical quantities 
like wave functions, energy levels, etc. The complete in- 
formation about these quantities can be extracted from 
distribution functions that can be calculated either nu- 
merically or analytically. Various methods of studying 
the statistics have been developed recently and this gave 
the possibility to achieve a good understanding of many 
disordered and ballistic chaotic systems (for a recent re- 
view see, e.g. 

Now it is well established that under very general 
conditions the statistics of chaotic systems is universal 
and well described by random matrix theory (RMT) or 
the zero-dimensional version of the supermatrix cr-model. 
For the distribution function of the amplitudes of wave 
functions one obtains within this approximation the fa- 
mous Porter- Thomas distribution [pi, which is simply a 
Gaussian distribution of the amplitudes. Then, fluctua- 
tions of the wave functions at different points are statis- 
tically independent. 

In this Letter, we present results from studying statis- 
tical properties of wave functions of two coupled chaotic 
systems. To be more specific, we consider two quantum 
chaotic billiards coupled in such a way that electrons (or 
electromagnetic waves) can penetrate from one billiard 
into the other. Studying such types of models is definitely 
interesting, because they can be relevant for double quan- 
tum dot systems recently fabricated on the basis of a 2D 
electron gas and investigated experimentally |7]] . One can 
also imagine that in the nearest future microwave exper- 
iments on double billiards will be performed. However, 
the importance of studying coupled chaotic systems is 
much more general and results obtained can be relevant 



e.g. for systems of weakly coupled complex atoms and 
molecules. To the best of our knowledge, neither numer- 
ical nor analytical results from a study of wave functions 
of coupled chaotic systems exist in the literature. 

Wave functions are the most interesting objects in cou- 
pled systems, because even at very small coupling they 
can drastically change if any two levels of the correspond- 
ing isolated systems are close to each other. Statistics of 
such hybridized wave functions can be quite unusual, but, 
at the same time, possess interesting universal features. 

The probability for the local density V \(j) a (x)\ 2 of a 
state with wave function (j> a (x) and energy e a as well 
as its spatial correlations can be described by the joint 
distribution function (JDF) f^ipiiP^), defined as: 



hipi,P2) 



where the brackets (. . .) denote disorder averaging, A = 
(vV)' 1 is the mean level spacing, V is the volume of 
the total system and v is the average density of states. 
The function /2(pi,P2) is the probability that the local 
densities at energy e are equal to p\ and pi at points X\ 
and X2 correspondingly. 

If pi and pi are not very large the statistics of the wave 
functions is universal and does not depend on whether 
the chaotic motion is due to disorder or a non-integrable 
shape of the billiard. The direct calculation using the su- 
persymmetry technique or RMT gives for single billiards 
in the unitary case: 



h(pi,P%) = exp(-pi) exp(-p 2 ), 



(2) 



which is just the Porter-Thomas distribution Proper 
formulae can be written for the other ensembles but we 
concentrate here on studying a chaotic system with bro- 
ken time reversal symmetry. This case is most simple 
for calculations but the interesting effects obtained be- 
low are general. The unitary ensemble is easily realized 
in quantum dots by applying a magnetic field. As con- 
cerns microwave cavities, time reversal invariance can be 
broken by a special preparation of the cavity walls || . 

The function /a (pi,P2), Eq. (^), characterizes correla- 
tions of the amplitudes of the electromagnetic waves in 
microwave experiments and can be measured directly || . 
At the same time, this function determines the statistics 
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of the conductance peak heights G max in quantum dots 
in the regime of Coulomb blockade. The conductance 
G max can be expressed through the tunnel rates of the 



contacts to external leads, I\ cx |</>(:Ei)| 2 flio| . p"l| |: 



G„ 



11 2 



t r : 



(3) 



Eq. (||) is valid for r, « T « A, where T is temper- 
ature and A is the mean level spacing in a single dot. 
Then we have for the distribution function P (g) of the 
dimensionless conductance g = T/(Te 2 ) G max : 



P(g) = J dpi d P2 6 (g 



P\Pl 
Pi +P2 



/ 2 (pi,p 2 ). (4) 



The relevance of Eq. (|^) for conductance fluctuations in 
single dots has been confirmed in experiments Un- 
fortunately, conductance fluctuations were not studied 
in an experiment on double quantum dots (?]]. In a re- 
cent theoretical work Jlj| the averaged resonance conduc- 
tance and its fluctuations have been considered for two 
coupled dots in the limit of weak coupling based on the 
assumption that wavefunctions of the corresponding iso- 
lated dots are Porter-Thomas distributed even near the 
point-contact. In contrast, we present a general micro- 
scopic derivation of the function f 2 (pi,p 2 ), Eq. ([!]), for 
two coupled quantum chaotic billiards for arbitrary cou- 
pling. We find novel universal correlations that do not 
exist in a single billiard. 



Assuming that the particles tunnel through a surface S 
between the two dots, the tunneling matrix element t a $ 
takes the form: 



tap — k 



id-1 



(8) 



In Eq. we integrate over the surface S, Vi is the vol- 
ume of dot i, to is the tunneling rate per unit surface and 
4> a (x) is the single particle wave function of state a at a 
point x of the tunnel contact in dot i. 

Starting with the Hamiltonian, Eqs. we perform 

calculations using the supersymmetry technique [Q . Fol- 
lowing this method the JDF can be written in a form 
given by an integral over supermatrices Q%, Q 2 [fblf : 

/a(pi,P2) = lim ~J o dtdQ 1 dQ 2 s(p 1 + ^Str(H, 6 Qi)) 
S (pa - 7(1 2A ^ Str(IUQ 2 )) C[Q t ] exp (-F [Q]) . (9) 



The free energy F = F 
for the single dots: 



Ft contains the free part F 



fo=-]E A" 1 Str(AgO, 



(10) 



and the term Ft describing the coupling between them: 

OL, 



Fn 



-Str(QiQa) 
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FIG. 1. Two quantum billiards are coupled along some 
extended tunneling surface S. 



The system we consider is represented in Fig. (|I|) . The 
Hamiltonian H can be written as: 

H = H 1 +H 2 +H T , (5) 

where Hi and H2 correspond to the isolated dots: 



a;i=l,2 



In Eq. (||), e l a are the single particle energies of dot 
and c + and c are creation and annihilation operators. 
We want to emphasize that the billiards are not identi- 
cal and can have arbitrary shapes and configurations of 
disorder. The coupling between the billiards is described 
by a tunneling Hamiltonian Ht that can be written as: 



l + „2 



C.C. 



(7) 



The function C[Q t ] in Eq. (§) is: 



0^} = Str ((n a/ - n r/ ) (D 1 Q 1 + D 2 Q 2 )) . 



A, = 
D l = 



A/Ai 



is 
and 



the 
A- 1 



mean level 



spacing in dot i, 
V2). The projec- 



tors liij select certain elements of the supermatrices Qi, 
A = diagjl, 1, 1, 1, —1, —1, —1, —1}, Str is the supertrace 
(for more details of the notations see 0). 

The parameter a cx \ d F T 1 S\t \ 2 (A^)" 1 is the di- 
mensionless coupling constant, with being the Fermi 
wavelength. In Ft we have kept only the term of low- 
est order in a. Other terms of order a n are smaller by 
a factor ~ (A^r 1 /^) and this corresponds to many- 
channel tunneling between the dots. Therefore, as long 
as we can consider not only small couplings 

but also a > 1. On the other hand, neglecting higher 
order terms for point contacts is correct only for aCl. 
Note, that in the limit a — > 00 fluctuations of (Qi — Q2) 
are suppressed and one comes to the model for a single 
dot with volume V = V1 + V2. The opposite limit, a 0, 
can only be performed after taking 7 — ► 0. 

To perform the integration we use the standard 
parametrization H . All manipulations are not very diffi- 
cult for the unitary ensemble. After a rotation of Q 2 only 
one set of Grassmann variables remains in the exponent, 
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which is convenient for integrating out the Grassmann 
variables. Due to the presence of the delta functions, the 
integration over the non-compact bosonic sector is sim- 
plified in the limit 7 — > and we get the final result valid 
for arbitrary values of the coupling a: 



hipi,P2) = Va exp {-Dipx - D 2 p 2 - P) 



(12) 
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Di{q k + a) 1/2 A l + B, l (q i + a) + C l (q l + aY 



In all the expressions i ^ k and both indices can take 
values 1 and 2. We have defined for compact notation 
qt = 2D k pi, (3 ■ 



A; 



2>D x D 2 m 



a)(q 2 
-(3D 2 



a) and: 
f D 2 ) n, 



B t = {Dim + ZD x D 2 n) (l + /T 1 ) f3~\ 
a = Bin (1 + 3/T 1 + 3/3- 2 ) /T 2 , 
sinh a 1 



cosh a 



sinh a. 



2a ' 2 

Eq. (p^), although being exact, is rather complicated 
and we discuss now its asymptotics. In the limit a — > 
corresponding to almost isolated quantum dots the JDF 
vanishes for all p\ , p 2 excluding the case when either pi 
or p 2 is zero, where it diverges. The divergencies are 
present in the terms proportional to /3 -4 in Ci and the 
terms proportional to f3~ 2 in B t . Then we find: 

f2(pi,P2) = D 2 2 8( Pl )e- D ^+Dl5(p 2 )e- D ^. (13) 

The particle is localized either in dot 1 , then p 2 = and 
Pi is distributed via Porter-Thomas, or vice versa. The 
factors Di appear, because the dimensionless densities pi 
have been normalized with respect to the volume of the 
total system and not of the corresponding single dot. 

In the opposite limit a — > 00 of strong coupling be- 
tween the dots we can easily obtain from Eq. ( |l2| ) the 
Porter-Thomas distribution, Eq. (|^), for the total sys- 
tem with volume V. 

So, in the limit of strongly coupled and isolated bil- 
liards we obtain the Porter-Thomas distribution. In both 
cases fluctuations of the amplitudes of the wave functions 
at different points are Gaussian and not correlated. But 
what happens if the coupling constant a is finite although 
small? Can one have for such a correlations between 
the wave functions at different points within one dot or, 
maybe, correlations between different dots? 

The answer to this question can be found taking in 
Eq. ( |l2| ) the limit a«l and pi,p 2 ^S> a, which gives: 

/a(pi,P2)» \^(D 1 D 2 f 2 exp \-(^+^f] (14) 



/2 



-1/2 



+ ( 2 + V U 1 U 2 ) ( u l + u 2 



3/2 



where ui = DiPi. We see, that the exponential cannot be 
factorized as in the case of the ordinary Porter-Thomas 



distribution, which demonstrates novel correlations that 
are absent in a single dot. This quite remarkable re- 
sult means that putting a weakly penetrable wall in a 
quantum billiard yields correlations of the wave functions 
across the wall, although without the wall they would be 
absent. Of course, making the wall less and less pen- 
etrable reduces correlations due to the prefactor ^fa in 
Eq. ( |l4| ) but this decay is slow. 

The form of the function f 2 (pi,p 2 ), Eq. (|i"4|), is sur- 
prisingly universal with only one sample specific param- 
eter a describing the tunneling. It is natural to assume 
that such universality follows from some general physi- 
cal principles. Indeed, the correlations between the wave 
functions of the chaotic billiards is a consequence of the 
quantum mechanical hybridization of states. This hy- 
bridization is most effective when the coupling between 
the billiards is weak and when any two levels of different 
billiards are close to each other. 



Let us show that Eq. (14) can be obtained in the limit 
a-Cl, pj.,2 3> 1 from this picture. The derivation is not 
as general and we are able to reproduce Eq. ([l4|) for point 
contacts only. Instead of Eq. (||) we consider a tunneling 
point contact at x t : 



(15) 



It 1 2 

The dimensionless coupling constant is now a oc • 
We assume that the main contribution comes from con- 
figurations where two levels of the dots 1 and 2 are close 
to each other and so we consider only two eigenstates 
with energy £ in dot 1 and energy £ 2 in dot 2. Including 
the tunneling results in a hybridization of the states. Fol- 
lowing p3[ we write the energies of the hybridized states 
in first order of degenerate perturbation theory: 



(16) 



where £+ = ±(£ J + £ 2 ) and £_ = i^ 1 - £ 2 ). The local 
densities | <^>± (ici .2 ) | 2 in dots 1,2 at the points x\ and x 2 
near the contact can now be expressed in terms of the 
eigenfunctions 1 (x\) and 4> 2 (x 2 ) of the isolated dots: 



l<M*i)l 2 
l<M^ 2 )| 2 



where A T = \t 12 \- 2 (£_ T J£ 



[1 + ^ 



1 



(17) 



tl2| 



In order to calculate the function f 2 (pi,p 2 ), the lo- 
cal densities |</>±(a;i ! 2)| 2 , Eqs. ( p~6| - |T7| ) , should be substi- 
tuted into Eq. ([!]) and the latter averaged over disorder. 
Instead of doing so, we replace the averaging over the 
disorder by an averaging over the level mismatch and 
the densities K = Vi\4> l (xi)\ 2 in all points. Assuming that 
the densities in the isolated dots are not correlated, be- 
ing Porter-Thomas distributed, and the level mismatch 
distribution being given by a function s: 
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(18) 



we can perform the averaging for the function f% (p\,P2) 
Now Eq. (0) becomes: 

f2ipx,n) = IE/" (ll d * e " 9i ) £l de - s(e - /A) 



As a natural application of the general formula, 
Eq. (12), let us calculate the distribution function for 



where A T = (tlqiq 2 ) ( £- =F \/ £- + *o<2i?2 



Evaluation of the integrals is straightforward. First 
we integrate over , 54 and then over £_ . In the limit 
V\tVi 3> a, a <C 1 the main contribution to the integral 
comes from £_ ~ io and we can replace the function 
s (£_/A) by the value s (0). This shows that the result 
is independent of the distribution of the mismatch 
The remaining integrals can be performed easily and we 



reproduce Eq. ( 14 ) up to a constant factor of order unity 
which cannot be found from such a consideration. 

Nevertheless, this simple evaluation cannot replace the 
explicit derivation of the general Eq. (|i"2|) with the su- 
persymmetry technique, because Eq. ( |12|) is valid for 
an arbitrary a and p\ t 2- Even in the limit a <C 1, 
Pi, 2 3> 1 the evaluation based on the two-level approx- 
imation was done for point contacts only. If we used 
instead of Eq. (|||) the more general Eq. (jg), we would 
have in the expression for |ii2| 2 combinations of the type 
J\ i (j> l *(x) <f> l {x') with \x — x'\ ~ Xp , However, the Porter- 
Thomas distribution is not applicable for describing wave 
functions at different points separated by atomic dis- 
tances. 

In contrast, Eq. (|l2|) does not depend on the structure 
of the tunnel contact and contains as the only parame- 
ter a. We have, therefore, identified universal statistics 
of hybridized levels that is independent of microscopic 
details and follows only from the existence of a coupling 
between two otherwise statistically independent systems. 



FIG. 2. Conductance peak distribution P(g) for different 
couplings a — 0.01 (dotted line), 0.1 (dashed line), and 1 
(solid line). 



the peaks of the conductance in a double dot structure 
in the regime of Coulomb blockade. Substituting Eq. jl^ ) 
into Eq. (^) we can evaluate immediately the distribution 
function P (g) , where g is the dimensionless conductance. 
The result of the computation is represented in Fig. (||). 
It should be mentioned that Eq. (^) and, hence, the plot 
in Fig. (||) are valid for a symmetric setup with both dots 
of the same volume and equal couplings to the leads (gen- 
eralizing to an asymmetric situation can easily be done) . 
For small couplings a <C 1 the statistical properties of 
the peak conductance are correctly described by the two 
level approximation In this region small values of 

the conductance peaks are most probable. For large a 
the distribution is almost flat. 

In conclusion, we derived the joint distribution func- 
tion of the amplitudes of wave functions for coupled 
chaotic systems in the unitary ensemble. We have dis- 
covered novel universal correlations that are absent in a 
single chaotic system and presented simple physical ar- 
guments explaining the origin of these correlations. The 
results obtained can be applicable for describing a large 
class of coupled chaotic systems. 
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